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Abstract 

Isolated electron spins in semiconductor nanostructures are promising qubit candidates for a 
solid state quantum computer, There have seen truly impressive experimental progresses in the 
study of single spins in the past two years. In this paper we analytically solve the Non-Markovian 
single electron spin dynamics due to inhomogeneous hyperfine couplings with surrounding nuclei 
in a quantum dot. We use the equation-of-motion method assisted with a large field expansion in 
a full quantum mechanical treatment. We recover the exact solution for fully polarized nuclei. By 
considering virtual nuclear spin flip-flops mediated by the electron, which generate fluctuations in 
the Overhauser field (the nuclear field) for the electron spin, we find that the decay amplitude of 
the transverse electron spin correlation function for partially polarized nuclear spin configurations 
is of the order unity instead of 0(1/N) (N being the number of nuclei in the dot) obtained 
in previous studies. We show that the complete amplitude decay can be understood with the 
spectrum broadening of the correlation function near the electron spin Rabi frequency induced by 
nuclear spin flip-flops. Our results show that a 90% nuclear polarization can enhance the electron 
spin T2 time by more than one order of magnitude in some parameter regime. In the long time 
limit, the envelope of the transverse electron spin correlation function has a non-exponential 1/t 2 
decay in the presence of both polarized and unpolarized nuclei. 

PACS numbers: 85.35.Be, 76.60.-k, 03.67.Lx, 
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I. INTRODUCTION 



The electron spin dynamics in semiconductor nanostructures is presently of particular 
interest both experiment ally^ 3 -* 4 - and theoretically because of the potential applications in 
spin quantum computation.— Among the many roadblocks to solid state quantum computing, 
maintaining the quantum coherence of the electron spin in a quantum dot is a crucial 
issue. One of the most important and relevant spin interactions for this confined electron is 
hyperfine coupling with the surrounding nuclei, or the so-called Fermi contact interaction. 
The role of hyperfine interaction could be both negative and positive depending on the actual 
experimental operations. On the one hand, nuclear spin bath could be a major decoherence 
channel for the electron spin that acts as a qubit.— ^ This problem is unavoidable in a GaAs 
quantum dot where all the isotopes of Ga and As has nuclear spin / = |, Although this 
is not a problem for silicon quantum dot where the nuclei of 29 Si (J = |, with natural 
abundance 4.68%) could be removed, so that the silicon is made up of only 28 Si or 30 Si, both 
having nuclear spin 1 = 0. On the other hand, electron- nuclei hyperfine interaction could 
also be utilized to control the nuclear spins, which as an ensemble could act coherently. For 
example it has recently been suggested that the ensemble of nuclear spins in a quantum dot 
may be used as a long-lived quantum memory for electron spin by transferring the electron 
spin dynamics to nuclear spin reservoir coherently.— 

The study of Non-Markovian electron spin dynamics in the presence of hyperfine inter- 
action is a complicated problem due to its quantum many-body (electron and nuclei) na- 
ture. The problem has been studied by many researchers previousfy . 9 i 10 i n i 12 i 13 i 14 i 15 An exact 
solution 9 - has been found for an electron interacting with the fully polarized nuclear reser- 
voir. Various types of approximations, perturbation theory^ or an effective Hamiltonian 
approach , 13 i 14 have been employed in dealing with this problem. Numerical studies , 10 i 12 i 13 
limited by the exponentially large Hilbert space, can only be applied to small systems, typ- 
ically of up to 20 spins. Here we analyze the problem starting from the exact Hamiltonian 
using the equation-of-motion (EOM) method by working in the Heisenberg picture. Usually 
the problem is studied using the master equations for the system density operator.— We de- 
velop a systematic large field expansion method to simplify the equations dramatically. This 
way we can address the full quantum mechanical problem analytically using well-controlled 
appr oximat ions . 
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A confined electron in a gated quantum dot generally has a nonuniform envelope wave- 
function and a relatively large radius, which leads to an inhomogeneous hyperfine coupling 
with the nuclear spins on the crystal lattice. The effective number N of nuclear spins in- 
volved is 10 4 — 10 6 depending on the actual size of the quantum dot. It is exactly this 
non-uniformity that causes both electron spin relaxation and dephasing. Physically this is 
because the electron acquires a different phase factor each time it interacts with a particular 
nuclear spin with random orientation. This effect is totally quantum mechanical, and is 
quite different from the dephasing effect caused by inhomogeneous broadening or ensemble 
average (T 2 *),— which can be removed by spin echo technique.— 1 ^ 

Electron mediated nuclear spin interactions have been studied for a long time, in both 
metals and in semiconductors . 19 > 20 Although at a finite magnetic field the direct electron- 
nuclear-spin flip-flop is highly unlikely due to the Zeeman energy mismatch, higher-order 
processes where electron spin does not flip are possible. For example, conduction-electron- 
mediated nuclear spin interaction [Ruderman-Kitterl-Kasuya-Yoshida (RKKY)] has been 
studied for a long time in both metals and semiconductors. Here we focus on the electron- 
mediated RKKY interaction between nuclear spins on single mediating electron spin. Specif- 
ically we focus on the backaction of this effective nuclear spin interaction on the mediating 
electron spin and consequent electron spin decoherence. Helped by a systematic large field 
expansion, we solve the full quantum mechanical problem analytically and reveal the crucial 
importance of the electron-mediated nuclear spin flip-flop processes in the decoherence of 
an electron spin. 

We note that virtual (high order) processes can be classified into two categories, one 
involving electron spin flip and the other without electron spin flip. In a high effective 
field (external plus nuclear field), these two types of processes have completely different 
consequences. In terms of energy the process with electron spin flip needs large extra energy 
to compensate the electron Zeeman energy, since the nuclear Zeeman energy is much smaller 
and can be neglected. Such processes (real processes for the electron), no matter whether 
they first order or higher orders, are highly unlikely for large fields. Not surprisingly a 
negligible decay amplitude of order O(jt) for both diagonal (relaxation) and non-diagonal 
(dephasing) components of the density matrix element of electron spin was found from these 
processes. In the second category, no extra energy is necessary since there is no electron 
spin flip between the initial and final states. It is found these kinds of processes indeed 
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have significant contributions to T2 using numerical simulation,— though for only 13 nuclear 
spins from an effective Hamiltonian. To illustrate the role of virtual flip-flopping of nuclear 
pairs, we shall compare the exact solution of fully polarized nuclei where no virtual nuclear 
spin-exchange is allowed with that of the case in which one of the nuclear spins is in the 
down state initially, assuming that all the nuclear spins point up in the fully polarized 
configuration. 

In general, exact solution cannot be found for a non-quadratic Hamiltonian. From the 
perspective of equations of motion, this means that there is an infinite number of equations 
when one attempts to evaluate the Green's function. Therefore, a decoupling method has 
to applied to cut off the series of equations and close the system. In the usual many-body 
treatment, this amounts to calculate some thermal averages using the spectral functions 
determined from the Green's functions self-consistently. However, such approximation in- 
volving ensemble average is not available in our present case, where we are interested in the 
real time dynamics including the coherent part and the decaying part of a single quantum 
many-body system. In our study we assume that the effective magnetic field (denoting the 
total of external and nuclear fields) is large so that the electron Zeeman energy is much 
larger than the single nuclear spin Zeeman energy, and the direct electron-nuclei flip-flop 
is not allowed energetically. This assumption greatly simplifies our discussions by selecting 
one group of equations, which give the major contribution to decoherence. Typically, an 
effective field slightly larger than the fully polarized nuclear field is enough for obtaining 
self-consistent result. Our study has experimental relevance since the full Overhauser field 
is about 4 Tesla in GaAs quantum dots. 

The electron spin relaxation induced by spin-orbit interaction is quite slow in QDs be- 
cause of the level discretization.— Nuclear magnetic dipolar interaction gives rise to spectral 
diffusion and also leads to dephasing at a time scale of 10 /xs . 14 i 22 i 23 The fluctuation of the 
nuclear field (the Overhauser field) is due to direct nuclear spin flip-flops. Generally the 
nuclear dipolar coupling is weaker than the hyperfine interaction in a semiconductor QD. 
Therefore, it would be interesting to compare the electron spin decoherence times caused by 
nuclear dipolar coupling and the virtual processes due to hyperfine interaction in the strong 
field limit. 

In this paper we present more detailed and extensive studies than those given in our 
previous paper.^ The remainder of the paper is organized as follows: in Section II, we at- 
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tempt to represent the electron spin decoherence, especially pure dephasing, with a properly 
constructed Green's function. We set up the equations of motion for calculating the cor- 
relation function. In Section III we calculate the Green's functions with several different 
nuclear polarization for the nuclear spin reservoir. The spectral functions which determine 
the decay behaviors of the electron spin correlation functions at large time are found. In 
Section IV and V we discuss the obtained results and conclude the paper. 

II. METHOD 

The hyperfine interaction between nuclei and electrons in the conduction band of a semi- 
conductor can be modeled by a simple Hamiltonian , 19 i 24 



Here the first term represents the Zeeman energy of the electron spin with uoq = g*/i#i?o; 
the second and last terms are the hyperfine interaction between the electron spin and the 
nuclear spins. Between the two, the second term represents the nuclear Overhauser field on 
the electron, while the last term describes the flip-flop of the electron and nuclear spins. The 
Zeeman energy of nuclear spins are neglected, because the Hamiltonian in Eq. (J1J) conserves 
the total spin angular momentum, and a single spin nuclear Zeeman energy is much less than 
that of the electron spin. Here A k is the hyperfine coupling constant at the fcth nucleus. It 
is proportional to the probability of the electron at position r and may be written as 



where ip(r) is the electron wavefunction (envelope wavefunction and Bloch wavefunction at 
the site of the nuclear spins). A is the total coupling strength between the electron and 
the nuclei: / A(r)dr = A. In a GaAs quantum dot, A m 90 \ieV . It is convenient to 
work with dimensionless units, and assume A/N = 1. This means that energy is measured 
in units of A/N and time is rescaled by N/A, with h — 1. A general form of A k will 
make the summation over fc quite complicated when looking for analytical results, especially 
for anisotropic three-dimensional quantum dots. A simple form, which is appropriate for 
two-dimensional gated quantum dots and has been used in previous calculations,— 



H = uj S z + £ A k F k S z + £ -±{1^8+ + liS~). 



(1) 



k k 




(2) 



is adopted in the following calculations. Here N is the effective number of nuclear spins in 
the dot. This coupling describes a two dimensional quantum dot with electron wavefunction 
decreases as a Gaussian as the distance from the center of the quantum dot increases. The 
index k labels the radial coordinate r 2 . A = J °° A^dk = N. In other words, A and N can 
be used interchangeably. 

The central quantity that we are going to calculate is the retarded Green's function of 
the localized spin state 

G ± (t) = -i9(t)(y \S~(t)S + (0)\y )- (4) 

Here 6(t) is the step function. \l/ is the initial wavefunction of the system (including electron 
and nuclear spins), 

|*o) = hH> + A)|ft}]|^n)- (5) 

We assume that the electron and nuclear spins are in a product state at t — 0, and are 
therefore not entangled. Furthermore, we assume that initially nuclear spins are in a product 
state \ip n ) = | If , If, • • • , i|, • • •), where II =f or |, i.e., nuclear spins could be either up or 
down at some particular site with a total net nuclear spin polarization P that will be defined 
later. Here we assume that the magnitude of nuclear spin is 1/2 although 1=3/2 for all the 
isotopes of GaAs. This assumption simplifies the algebras in the following study though we 
do not anticipate this simplification to cause any qualitative difference in the properties of 
Gj_(t). A mixed nuclear spin state could be expressed in terms of a linear combination of 
the product states. Thus, the results of our calculations of product states can be used to 
study more general initial nuclear states directly. 

An important quantity that represents the quantum coherence of the electron spin is the 
quantum mechanical expectation values of transverse electron spin operators, (^ol^'WI^o) 
or similarly (\E f |S' + (£)| 1 I f o)- The Green's function defined in Eq. (@J) is equivalent to 
(v&ol-S^WI^o) in the high field limit. Substituting Eq. © into the expression (^ |5 ,_ (t)|^ ) ) 
we find 

(tfol^WW = a*(3 (^ij n \e im S-e- iHt \ fafa) + higher order terms. (6) 

Among higher order contributions, one or more electron spin flips can lead to non-zero 
results. However, as we have argued in the introduction that these terms are of the order 
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0(1/N) or 0(1/N n ) for polarized nuclei in the absence of external an external magnetic field. 
The leading order contribution in Eq. (jUJ) on the contrary does not involve any electron spin 
flip during the evolution. In fact this term gives rise to pure dephasing (as the electron spin 
does not flip) while the higher order terms involving electron spin flip lead to decoherence 
induced by relaxation. The leading order term in Eq. (jUJ) is exactly the Green's function 
defined in Eq.(@J) up to a proportional constant with the assumption that |\l/o) = I ^,ipn)- 
Another way to understand the pure dephasing term in Eq. © is to consider the phase 
evolution of the electron spin in the Schodinger picture, 

G ± (t) = -t9(t)(^F kir --J^e im / h S~e~ UH / h S + \^F kl ,---,F kN ) 

= -i6{t) {(4; I z kv ■ ■ ■ , I z kN \e im > % } S~ {e im ' n \ - - - , F kN )} . (7) 

The term in the first curly brackets represents the evolution of the electron spin-down state 
in the presence of the hyperfine interaction, while the term in the second curly brackets 
represents the evolution of the electron spin- up state in the same environment. If no electron 
spin flip occurs, any decay in the calculated average will be due solely to dephasing between 
the electron spin-up and spin-down states. Obviously, electron spin flip will also cause decay 
of the correlation function. Therefore, G±(t) contains the complete decoherence information 
for the electron spin in consideration. 

Without loss of generality, we assume |^ ) — I •l^V'n) m the following discussions. The 
equation of motion can be set up in the Heisenberg representation 

i^G ± (t) = S(t) - i9{t){V \[S-(t), H}; S + (0)|v|> >. (8) 

The solution of this equation can be obtained in the energy space by performing the Fourier 
transformation, after which we have an equation in the form of the standard equation of 
motion for the correlation function of two arbitrary operators 

u((A; B)) u = (tfo|i£|W >| t =o + «[A H]- B)) u , (9) 

where expectation value in the first term on the right hand side is calculated at the initial 
time, and ((A;B)) U = J(-i)6(t)(ty \A(t)B(0)\ty )e iuJt dt. 

Although the hyperfine Hamiltonian [Eq. ((TJ] looks quite simple, the general solution has 
not been found. Instead one could obtain approximate solutions under various conditions. 

7 



Exact solutions can be found in two simple cases: namely the system with fully polarized 
initial nuclear spin configuration, and the Hamiltonian with uniform coupling constants. 

It is straightforward to derive the following equations of G±(u) by calculating the spin 
commutators repeatedly, 



k k 

( u 2 -^)((I k -,S + )) ul = -^ + ^)G L {u J )+u J A k ^n k S~ ] S + )) w + \Y.A k A kl V kk ,{u) 
4 1 1 1 k'(k) 



(u?-£)mS*-,S + )) u = ^((n k S--S + )%-^{u + ^)G L {u) + uY.^V kk ,{u), (10) 

fc'(fe) 

with 

Vkk'{uj) = ((I k S + I«; - ((I k I+ t S-; S+)U (11) 

where n k = | — I k and Q = u + y. Here Q is the effective magnetic field for fully polarized 
nuclear spin reservoir. J2 k '(k) represents the summation over k! for k! ^ k. After rearranging 
the equations, we find 



[w - fi - Eo(w)] G L {u) = 1-EM^;S%-^7T^M^M (12) 

k 1 k^k' \ u ~ A klV 



with the self-energy 



We will use Vt = cu + Yl k A k I k to represent the effective magnetic field experienced by the 
electron in the following discussion. The meanings of the various correlations are obvious. 
For example, ((n k S~; S + )) UJ represents the Fourier transform of the time-correlation function 
for flipping the electron spin to the down state while the kth nuclear spin is also in the down 
state {n k = 1). Both ((I k ; S + )) UJ and ((I k ~S z ; S + )) UJ describe the direct electron-nuclear-spin 
flip-flop. The only difference is that the latter one depends on the electron spin state. 1/Q 
will be taken as an important expansion parameter in next section where we look for solutions 
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of the Eqs. (fTUj). Q ~ N corresponds to polarized nuclei or external magnetic field with a 
few Tesla's. Q ^> 1 will be a good starting point for simplifying the EOMs. ((rikS~; S + )) UJ 
and Vfcfc'(cu) are higher order correlation functions. We should mention that the Eqs. (|1(J|) 
are exact. Approximations will be applied to obtain the general solution of arbitrary nuclear 
polarization when Q ~ N. We will repeatedly use two different approximations, large TV 
expansion and large field expansion. These two terms have different meanings. The first 
one is just a mathematical treatment, where we simply neglect Ak when Ak and Q appear 
together since Ak ~ 1. This approximation simplifies our EOMs significantly. The large field 
expansion, which will only be used for solving the high frequency solution of the partially 
polarized and unpolarized nuclei, is a physical approximation, where we select a group of 
EOMs that describe the major contributions to the self-energies. 

III. SOLUTIONS OF THE GREEN'S FUNCTION 

In this section we obtain the solutions G±(u) of Eq. (fT0|) presented in Section II under 
various conditions. The real time dynamics of the Green's function, Gj_(t) can be found 
using the spectral function p(uj), 

/oo 
p(uj)e~ iuJt duj, (14) 
-oo 

where the spectral function p{u) is obtained by analytical continuation (u — > uj + i0 + ) 
using the retarded Green's function G±(u)^ 

P {uj) = --lmG L (u + iO + ). (15) 

7T 

The imaginary part of the retarded Green's function contains contributions from poles and 
branch cuts.-^ More specifically, the poles will determine the renormalized energy or the 
oscillation frequencies of the electron spin in our case; while the branch cut describes the 
decay of the electron spin state. Our focus in this section is to calculate these quantities for 
various nuclear spin polarization, with or without external magnetic field. 

A. Fully polarized nuclei reservoir 

The solution to the problem of an electron spin interacting with a fully polarized nuclear 
spin reservoir is quite straightforward. One crucial argument we will use to simplify the 
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EOMs is the conservation of angular momentum at any time during the evolution. For 
initially fully polarized nuclei (all nuclear spins in the up state) and initial wave function 
= I 4) "0n) j this means that there is only one spin (either the electron spin or one of the 
nuclear spins) in the down state. This restriction eliminates all the higher-order correlation 
functions present in Eq. (fT2|) . For example, ((n^S 1 "; S + )) u = because when the electron 
spin is flipped down, = 1/2 — 1% must be zero for fully polarized nuclei. Vkk'{uo) vanishes 
for the same reason. Therefore Eq. ()12j) can be reduced to 

G ^> = u - a ,'-E, H ' < 16 > 

where Q = ujq + N/2 for fully polarized nuclei. Assuming that the hyperfme coupling 
constant takes the form in Eq. Q, the real and imaginary part of the self-energy are 



1 r°° Ai N 1 

ReEnM = - / k —dk = Nuln 1 

2to 



(17) 



and 

7T f°° / Au\ 

ImE (cj) = / AlSluj )dk = -Nttu, for < u < -. (18) 

4 jo V 2 / 2 

In both calculations we convert the summation over k into an integration, and use the 
formula 

in finding the imaginary part. The imaginary part of the self energy is nonzero only when 
| > uj > 0, this is due to the constraint of the delta function in Eq. (fT%|). The spectral 
function can be calculated with Eq. (|TK|) 

Z p 5(u — f2 — ReSo) uj ~ N 

P(f*)=< ! 1 (20) 



Nuj r±s.,j. +ln i 1 _j_|i T l^ 



< w < |, 



where the renormalization factor Z p is^ 



9a; 



(21) 



The pole u p is the solution of equation uo — Q — ReS (c^) = 0. When ujq = (zero external 
field), u p « JV/2 + 1/4 and 2 P ^ 1 - 1/JV. It is clear from Eq. (JH that there are two 
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FIG. 1: To obtain the asymptotic behavior of G±(t)&t large time, the Fourier integral in Eq. (|14|) 
is converted into a Laplace integral by deforming the original contour C into C1 + C2 + C3 as shown 
above and then allowing uj — > 00. 

contributions to the spectral function p(u), namely a pole near Q and a branch cut in the 
range of [0, -]. For a schematic, see left panel of Fig. (J2J). In the fully polarized case, the 
weight of the pole is of order 0(1), while only a small fraction of the spectral function is a 
continuous spectrum. Since only the branch cuts or continuous spectrum contribute to the 
decoherence, electron spin in a fully polarized nuclear environment has an infinite T 2 time. 

In Appendix [X] we show that this exact solution can also be found with the Jordan- 
Wigner representation by reducing the original spin Hamiltonian into the non-interacting 
Anderson impurity models The physical meaning from the exact solution of the Anderson 
model is that part of the time the electron spin (localized state) is in the spin-down state, 
and part of the time one of nuclear spins (continuous state) occupies the down state. The 
true eigenstates (resonant scattering) are the linear combinations of these two states. 

In order to obtain real time dynamics we need to perform an inverse Fourier transform. 
This is straightforward for the delta function part of p(w), the contribution of the pole. 
However the branch cut integration is much more complicated, and a closed form cannot be 
found. To investigate the long time decay of the Green's function, the method of the steepest 
descent^ is used. The approximate calculation of the integration is justified by carefully 
deforming the integration contour in the complex u plane so that Laplace's method can be 
used. In Laplace's method the integration of an exponential function is approximated by 
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the integral of this function in the neighborhood of the global maximum of the integrand. 
In order to use the method of steepest descent to calculate the frequency integral in Eq. 
(|14|). we deform the original integration contour C into C±, C2, and C3 as shown in Fig. 
n The integration along C3 at infinity in the lower half plane is exponentially small in the 
large time limit. This is exactly the spirit of steepest descent, as one can see the term e~ luJt 
is negligible if u is analytically extended to the complex plane through u = —is and s is 
allowed to approach positive infinity. Substituting the spectral function into Eq. (jT4"j) we 
obtain 



G±(t) 











N Jo 



1 



-iu)t 



duo 



£ + ln 



^ 2w 



(22) 



where we have neglected the factor —i6(t) in Eq. (jl4j) and assumed that c^o = 0. On the 
contour C\ where u = —is, the integral is determined by the integration interval near s = + 
as t — > +00 because of the term e~ st . By performing a Taylor expansion around s = + 
we obtain p±(—is) = — -j—. The asymptotic form of the integral for the inverse Fourier 
transform when t — ► 00 is 



G x (t) 



1 

+ N 



+ O 



h ( ln2 - ln7 



71 

i— 
2 



1 



(mt); 



N (\ntf 



(23) 



where 7 = — / °° lnse s ds is the Euler constant. It can be shown that the integral along C3 



has a smaller contribution that is proportional to 



1 



Ntlnt ' 



The amplitude of the decay part 



is only of the order O (h\ an d the leading order behavior at long time is l/\nt. We have 
written the oscillatory part as [1 + 0(-^)]e l '^ t , where we have neglected the 1/N correction 
of the coefficient. The results thus agree with those found previously.™ 



B. One-spin-flipped nuclear reservoir 

The exact solution for the fully polarized nuclear reservoir is instructive but of little 
practical relevance in current experimental situations, where maximum achievable nuclear 
polarization in semiconductor quantum dots is in the order of 50%. A pertinent question 
is whether there is any qualitative difference between partially polarized and non-polarized 
nuclear reservoirs in the spectral function, spin correlation, and its decoherence To address 
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this question, we consider another simple case in which there is only one flipped nuclear spin 
in the otherwise fully polarized initial state, i.e., |W ) = | fi, T2, ■ • • , Ik, ■ • • , Tat), where k 
is the index of the initially flipped nuclear spin. 

The EOMs are quite complicated even for this simple case. In Appendix [B] we give the 
complete EOMs, including all the higher-order correlation functions. Together with Eq. 
(I12J1 . these six equations form a closed set. The even higher-order correlation functions 
vanish because of conservation of angular momentum, as we have discussed in the case of 
the fully polarized nuclear reservoir. We point out that these equations are exact. They are 
obtained from the iterative EOMs by calculating the commutators. It is a formidable task 
to find the exact solution of these equations, thus we first simplify the equations as much as 
we can. It turns out to be convenient to discuss the low energy (\u\ < |) and high energy 
{u ~ Q) solutions separately. From the exact solution of the fully polarized case, we already 
know that the spectral function behaves quite differently in these two regimes. Specifically, 
the low energy part corresponds to a continuous spectrum, while the high energy part is just 
a delta function that gives rise to the undamped coherent oscillation in the electron spin. In 
the following we shall discuss both the low energy and high energy solutions of one-flipped 
nuclear spin reservoir with 1/N expansion. Since N is very large, this approximation should 
make little difference from the exact solutions. 

Low energy solution. To obtain the low energy solution, we start from Eq. ()12j) . By ex- 
amining the right hand side of Eq. (fT2j) . it is easy to see that one needs to find ((rikS~] S + )) UJ 
and Vkk'(u)) [Recall that Vkk'(u) = ((I^S + Ij7,; S + )) LU — ((I^I^,S~; S + )) LU ). To simplify the 
equations, we use the fact that 10 <C Q in the low energy regime and <C Q. As we point 
out in the introduction, we assume the existence of a large effective field so that ~ N. 
Therefore we can use the approximation, u ± f2 =1= j4fc "^ j4fc/ « ±fl In other words, we can 
neglect u and A/, whenever they appear in a sum with Q. This approximation helps find 
the leading order solutions. From Eq. (Bl), we obtain ((nkS~; S + )) UJ , 

((n k S-; S + ))„ = - * , > *)o + 9 " ? T] E -V^WM + (4) • (24) 

Here the self-energy term T, (uj) is of the order O(N) because it involves a summa- 
tion over nuclear sites [also see Eqs. (|18|) for the fully polarized nuclear reservoir]. By 
studying Eq. (|10|) and the equations in Appendix El one can draw the conclusion that 
Vkk'(uj) ~ ((rifciS' - ; 5' + }} a; . This is because Vkk'{w) ~ {{1^ S~; S + )) w /Vt, which is a result 
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of Eq. (B4). By putting this relation into Eq. (|10|) one could immediately see that 
Vfcfc'(u;) ~ ((rikS~; S + ))uj. With these results we can determine the order of magnitude 
of the summation J2k Ak((nkS~; S + )) w using Eq. (|24|) 



1 



^2(n k ) oc O 



(25) 



k * u 1 ^oWv fc 

since J2k( n k)o — 1- [ Recall that there is only one nuclear spin in the down state (n^ = 1) 
initially, and all other nuclear spins point up (n k = 0) ]. Similarly one can easily find that 
both the third and the fourth terms on the right hand side of Eq. ()12|) are of the order 
O (fi)- Therefore we obtain the leading order low energy solution of G±(uj) from Eq. (fT2*|) 



G] M 



O 



(26) 



for uj ~ 0(1). This solution is exactly the leading order approximation of the fully polarized 
nuclei in the low energy limit. Clearly there are only higher-order differences for the two 
cases in this regime. The leading order spectral functions also have the same behavior. 

High energy solution. The approximate solution in the high energy limit is much more 
difficult to obtain compared to the low energy limit. In this limit uj ~ f2, thus A k <C u, fl. 
Generally we can drop A k in the expressions such as uj ± 4r and ± 4^-. However, we must 
be careful when dealing with terms like uj — Vt ± 4^-. Since both uj and f2 are of order N, 
the difference between them could be 0(1). Thus we can no longer neglect the A k terms in 
these expressions. With Q ^> Aj,, we find Eq. (IB2|) can be simplified to 



({n k S~-S + )) L 



1 



UJ 



(n k )i 



4Q 



(fk + 9k) 



(27) 



Vt — Eo(w) 

where fk(uj) = J2k" Ak"Vk»k(w) and gkiyj) = J2k" Ak»Vkk"(uj)- Substituting Eq. (}2*Tj) into Eq. 
(fT2*|) . we obtain 



G] (uj) 



uj 



Q — Tmiuj) 



A k A k ,V w (uj) 



k^k' 



(28) 



Here Vkk'{uj) contains two terms, one of which is negligible compared to the other in the 
high energy limit. It can be shown that ((1^ S + I^; S ,+ )) u; ~ O \jA ((I^I^,S~; This 
relationship is to be expected as ((I k ~S + I k ~r 1 5 ,+ )) u; corresponds to the electron being flipped 
during the evolution, which is a probability event O (j^ . Similar simplifications will be 
seeked when we search for solutions of partially polarized and unpolarized nuclei. 
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In order to evaluate V k y(u), we first simplify Eqs. (jB3|) and (|B5|) using the large N 
expansion method we have outlined above, it can be shown that the two equations can be 
approximated with 

((I k ny, = -^{(n k >S-; S + ))„ + ^V w (u), (29) 

and 

((/ fc -/+J fc -,; = ^ W«) + ^V w (u). (30) 

Using Eqs. (jB4|) . (|2*9j) . and ffH0|) . we can derive an implicit equation for V k y 

A k + Ay\ A k Ay (n k ) A k A k > (ny) 



I oj H I VW(^) 



4fi - - S + A fc 4fi - - S + Ay 
+ ^[Afc^H+Afc// fc (u;)], (31) 



where Co = uo — Vl — ^ and A 2 = X^^i- Using the wavefunction in Eq. (J2J), it is easy 
to show that A 2 = f . Multiplying both sides of Eq. (J3TJ) by A fe // + (A fc + A fc >)/2] and 
summing over k', we find the equation for f k {oj) takes the form 



2 



, w afc-4fc(^fc)o ^fc^s 

(12 - a k )fk =~ 



^ fe"(fc) ^ + - 



where 



<■*(*> -jE^rfc- (33) 

fe" ^ H 2 

The equation for g k is obtained by interchanging f k and g k in Eq. (|32j). Through the 
symmetry of the equations for f k and g k , we find f k = g k in the leading order large Q 
expansion. From Eq. (J32J) it can be seen that J2kA k f k ~ 0(1) near the pole of l/(u;— f2— S ). 

Equation (|3*2*|) can be transformed to an integral equation by converting the summations 
into integrals, though this integral equation is still difficult to solve. For our present purpose, 
Eq. and Eq. (J32j) are enough to further our discussions. First, by solving the equation 
uo — Q — S (co>) = 0, we obtain the same pole as the fully polarized case. The interesting 
point here is that the pole now has a reduced weight, because the residue is Z p ps 1 — 
n ~~ ^.^k-^-kfk- Since J2kA k f k ~ 0(1) when to — Q — S = 0, the difference from the 
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FIG. 2: Schematic of the spectral function for fully polarized nuclei (lefty panel) and one-spin 
flipped nuclei (right panel). The shaded area in the right figure denotes the contribution of the 
branch cut by considering virtual nuclear spin flip-flop which is not possible in the case of fully 
polarized nuclei. Also the delta function represented by the vertical line has a slightly reduced 
weight compared to that in the fully polarized case. 



fully polarized nuclei is very small (Z p ~ 1 — 1/N in the exact solution). On the other 
hand, the function /fc(a>) has non- vanishing imaginary part when the frequency is in the 
neighborhood of Q. In other words, the spectral function has a branch cut near Q. This 
contribution to the spectral function is only of the order O f^J because of the prefactor 
^ in Eq. (|28|h Combining these discussions with the picture we have established in the 
low frequency region, we can summarize the properties of the spectral function in the right 
panel of Fig. 121 In the left panel, we draw the schematic of the spectral function of the 
exact solution in the fully polarized case for comparison. It is quite easy to understand 
where the new feature comes from. In the fully polarized case, the only possible process is 
the electron-nuclei flip-flop (see the upper figure in Fig. EJ). The electron flipped state is 
a high energy state, where energy hfl (neglecting the nuclear Zeeman energy) is needed to 
compensate for the Zeeman energy mismatch. For a one-spin-flipped nuclei reservoir, virtual 
nuclear spin flip-flop mediated by the electron is possible (see the lower figure in Fig. OJ). In 
the initial and final states, the electron spin state remains unchanged, but two nuclear spins 
exchange their states. The virtual (intermediate) state is a high energy state, and the cross 
section of this process is proportional to I/O. 



1(3 



C. Partially polarized and unpolarized nuclear reservoir 



The solutions of both the fully polarized nuclear reservoir and one spin flipped nuclear 
reservoir, while instructive, have little experimental relevance. It is important to find the 
solution for a general nuclear spin configuration when studying the electron spin decoherence. 
To proceed, we first define the effective nuclear polarization P by 

P = (34) 

where (iVj) represents the number of nuclear spins with up (down) spin. We further 
assume that both iVj and iVj_ are large, i.e., TVf ~ 7Vj_. We can then convert the summation 
over both the up and down nuclear spins into integrals. We expect the solution to be similar 
to those we have already found in the limit of highly polarized nuclei when N± <C iVj or 
Ni > iV T . 

We treat the nuclear field (or Overhauser field) J2k ^k^k with the adiabatic approximation. 
With a large effective magnetic field Q, the electron spin precesses with frequency S. On the 
other hand, the nuclear field fluctuates around its average value by a small amount in a much 
slower time scale, since there are many nuclear spins (N ^> 1) interacting with the electron 
spin simultaneously. Each of the nuclear spins has only a 1/N probability to exchange spin 
with the electron. Under this approximation we can neglect the small time-dependent change 
of the nuclear field Ylk ^k^k > an d factor it out from the time-dependent correlation function 
((J2k AkIk(t)S~ (t); S + )) . Meanwhile, the nuclear spin raising and lowering operators are 
treated fully quantum mechanically. This approximation reduces the higher order correlation 
function ((rikS~; S + )) UJ to G±(ou). The term YJ fc A^P^S Z in the Hamiltonian remains as an 
operator for calculating commutators. When we take the nuclear operator in the correlation 
functions as a c-number, we always make sure that the difference from the real quantity only 
gives rise to higher order contributions. 

Similar to the highly polarized situations discussed previously, we study the solution in 
the low energy and high energy limit separately. Recall that in the case of polarized nuclear 
reservoir with one flipped spin, the weight of the delta function in the spectral function is 
reduced by a small amount, and there appears a new continuous contribution to the spectral 
function near Q because of the nuclear spin flip-flop mediated by the electron. Since there is 
only one nuclear spin in the down state initially, scattering channels for nuclear flip-flops are 
limited (see Fig. EJ). Therefore, decoherence effect of the continuous contribution is small, 

17 



(a) 




fully polarized nuclear spin reservoir 




intermediate high energy state 



initial state 



final state 



FIG. 3: Schematic of possible spin exchange processes in case of fully polarized nuclei (upper panel) 
and one-spin- flipped nuclei (lower panel). For the fully polarized nuclei, the only possible process 
is the direct electron- nuclei flip-flop, which requires a large assisted energy to compensate in a high 
magnetic field. Second-order virtual processes are also accessible for the one-spin-flipped nuclei. 
In this case, the high energy state is an intermediate state, while the electron spin state is not 
changed in the initial and final states, and two nuclear spins exchange information effectively. 



of order O . With more nuclei in the down spin state, we expect that the continuous 
contribution will become increasingly important, and be of the order unity eventually when 
iVj_ ~ Np We will verify this conjecture in the following discussions. 

Low energy solution. Helped by the adiabatic approximation for the nuclear field, we can 
write the first equation in Eqs. in the following form, 



Thus the leading order contribution to G±(u) is O(^). We can also write the equations for 





(35) 



<<4-;S + )> w and((4-^;S + >> w as 



A k ({ItS*;S + )) u -h k G ± (u), 



(36) 



cu((I k S*;S + )) 



(37) 
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with hk = Afc(If)o an d V = S + J2k Aklk ~ Sfc Akl~£ S~ . By introducing hk, we are taking a 
mean field average for the Overhauser field produced by this individual nucleus. Substituting 
Eqs. (JHtj|) and (JBTj) into Eq. (J3T)|) . we obtain 



(38) 



Fortuitously, further calculation of {{I^V; S + )) LU does not lead to any other higher or- 
der correlation functions in the low u limit. This is illustrated by evaluating both 
((I k ~S + I~; S + ))u and ((I^I + S~; S' + )) w with the EOM by a 1/N expansion. Here we have 
defined 1^ = Y.k^t^k as collective nuclear spin operators. The index is dropped when 
n — 1. Neglecting all higher order terms, we find 



- Of)o«4- Q - S*) ; + (39) 

where we have used the operator equality S + S~ = | — S z and the fact that 

((Ij:S + I~;S + }} U} ~ O [{{I^S + I 2 ;S + }} LU ]. Combining Eq. with flUD, we find a sin- 
gle equation for ((I k ~S z ; S + )) UJ and G±(lj), 



P4, 

+ ^E^((4^;S + >>.. (40) 

In this equation, we have used the mean field average and, replaced the nuclear spin expecta- 
tion value of the z component (I^.)o by the average nuclear polarization 4j. The summation 
Y^kAl((IkS z ; S + )) UJ in Eq. pUj) can be found by multiplying both sides of the equation by 

A 2 

lu 2 -A 2 /4 an d then performing the summation over k. We then obtain 

- + %"° a ' )G± - (41) 

with 

.4 



E— "V (42) 



fe - 
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FIG. 4: To obtain the asymptotic behavior of Gj_(t) for partially polarized nuclei at long time 
limit, the Fourier integral in Eq. (|14j) is converted into a Laplace integral by deforming the original 
contour C+ and C_ into C±, C2, C3, and C4 as shown above and then allowing uj — > 00. 



a n ~ N since it is a summation over all the nuclear spins. Plugging Eq. (|41j) back into 
Eq. (J4"U|) . we obtain an expression for ^2 k A k {{I k S z ] S + )) u , which can then be used to find 
G±(uj) using Eq. (J35j) . After all these evaluations, we obtain the low energy solution for 

n i \ m - 2uja * + P(J 3 fAO\ 

G±{LJ) - -\8Q + Pa 3 y-P^ 2 (2N + a,y (43) 

The final solution is quite complicated. One important result is that G±(u) ~ O(jf), just 

like in the highly polarized cases. In this solution the higher-order correlation function has 

no effect since the term {{I k ~I + I~; S^))^, which involves nuclear spin flip-flop, does not 

appear in Eq. (|4()|) . The summation a n (uj) can be calculated with the same technique as 

we have used for calculating T,q(u) with analytical continuation (u — > uj + i0 + ) and the 

conversion of summation to integral (J2k ~ > Jo°° ^ ne rea ^ anc ^ imaginary parts of a n are 

evaluated and given in Appendix The real parts of a n have two branch cuts, one from 

to |; the other one from — | to 0. The imaginary part of G±(u) is non-zero only in these 

regions. 

We first consider the case of zero external magnetic field, when f2 = To study the 
long time behavior we again deform the original integration contours C + and C_ into Ci, C2, 
C3, and C4 in Fig. 0] The integrals at minus infinity in the lower half plane are neglected. 
Using the results of a n {uj) given in the Appendix we find the spectral function behaves 
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as 



x i(P±l) 1 

on contour C2 (u = 0~ — is) and C3 (a; = + — is) as s approaches + . Calculating the 
Laplace transform, we then obtain the long time asymptotic form of time dynamics of the 
Green's function in the low energy limit, 

G ^ K TnL> ^ °°' (45) 

The long time asymptotic form 1/lnt is the same as that of the fully polarized case, and 
we do have the right numerical limit as P — > 1. Consider now the situation of the strong 
Zeeman field limit where \ujq\ ^> N. The Green's function is proportional to ^ when w<iV. 
The asymptotic behavior of G±(t) will be 7, and the decaying amplitude is proportional to 
^. Again this result is the same as the fully polarized nuclei case.^i 

The properties of the low energy solution to the spin correlation function studied in this 
subsection are quite similar to those that have been investigated in previous studies.—^ The 
basic result is that the decay amplitude is of the order O (j^ in the perturbative limit when 
the effective magnetic field is comparable to the number of nuclear spins in the QD. However, 
only the direct electron-nuclei flip-flop has been considered in these studies. Therefore it is 
quite reasonable that the contribution is small since the Zeeman energy mismatch is large. 

High energy solution. In the following discussion we search for the solution to the spin 
correlation function in the high frequency region. This solution will be of the most experi- 
mental relevance. In the case of a highly polarized nuclear reservoir the major contribution 
to the spectral function in the large u> limit is a delta function, that does not lead to any 
decoherence in the electron spin. When the effect of the electron mediated nuclear spin ex- 
change is included, decoherence with a small amplitude arises. The question is whether this 
effect will become large enough so that there is a measurable contribution to electron spin 
dephasing when initially the two nuclear spin species (up and down) have approximately 
equal population. We address this question in the following discussion of the high energy 
solution for partially polarized and unpolarized nuclei. 

By using the adiabatic approximation and considering that uo = Q + 0(1) and uo 3> Af., 
Eq. fll2jl can be straightforwardly simplified to 

u - - £■) G ± {u>) = 1 - 4 E A k A k ,((Ij;l+,S-; (46) 



AO, J 2f2 k ^ k , 
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where A 2 = Y,k^l- Notice that 'S + I#; < ( (1^1^ 3~ ; because 

({I^I^iS~; S + )) LlJ is proportional to CJ _ n+ ( A 1 fc _ Afe/ - ) / 2 , while ((I^S + I^,; S + ))u is proportional 
to ^tq- Mathematically, the difference originates from the commutators of with S 12 , 
i.e., [S^S*] = Since = Q + 0(1), it is clear that a; _ n+(A 1 fc _ Afc , )/2 ~ 0(1) while 

ZJ+Ti ~ Physically this difference is also easy to understand. In ((I^S + Ij7,; S + )) u 

an electron spin needs to be flipped down, while in ((Ij7l£,S~; 5 ,+ )) a; , the electron spin is 
flipped up and then flipped down, and two nuclear spins at k and kl flip-flop with each other. 
Therefore the process in {{I^l£,S~; S + )) UJ does not require a real electron spin flip and is a 
nuclear flip-flop process mediated by the electron. From Eq. I|46j) it is clear that this is the 
only higher-order correlation function that contribute to Gj_(uj). 

The calculation of the one-pair correlation function ((I k l£,S~; S + )) UJ is non-trivial, since 
its computation involves two-pair correlation function ((I^I^I^I^S"; S + )) w , which in turn 
depends on even higher order correlation functions. We thus need a cut-off scheme to close 
the EOM. This is not a problem for the low energy solution where the terms of the one- 
pair correlation functions cancel each other so that the transverse electron spin Green's 
function only depends on the electron-nuclei spin-flip correlation functions S^))^ and 
((IfrS z ; S + )) UJ . In the high energy regime, it turns out that the self-energy function can be 
expanded in powers of Therefore, as long as the effective field fl is sufficiently large 
compared to N (meaning that the applied field is larger than the Overhauser field from a 
fully polarized nuclear spin reservoir), the approximate solution can be obtained with the 
first few terms in the expansion. 

We start with the exact EOM for ({I^I^S~; S + )) w , 

A h ,( 1 - + Il){(I^S z ] S + ))„- ((J fe -J+ £ A k nI w ,S z ;S + )) w . (47) 
J k"(k,k') 



The leading order (in terms of ^) EOM for ((I k Iy J2k"(k,k') Ak"Ik»S z ; S + )) UJ takes the form 

4 



E A k ,Jv f S z ;S + ))„ = -^((I^I+S--,S + )) b 



k"(k,k') 

((/+ E Ii>S-;S+)) m (48) 

k"(k,k') 



where two-pair correlation functions have been neglected. In leading order of 1/ft, the EOMs 
of the two lowest-order correlation functions of electron-nuclei flip-flopping in Eq. (jlUj) are 
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simplified to 



((Tj:- 1 ^)) u = ~A k G ± (u) 



and 



Ak 



G 



2Q 



((I k Y,Ik'S-;S + )l 



(49) 



(50) 



fc'(fc) 



Combining Eq. (J4T|) with Eq. (J5U|) we find the following equation for ((I k I k ,S ; S + )) u 



to 



A k -A k 



1-P 2 

8fi 



A k AuG i fwl 



+ 7j£« J * E ^^";^»« + ^<W E ^i^5'-;^» w . (51) 

4iZ k"(k,k') 4ii fc"(fc,fc') 



with u) = u — Q — which represents the deviation of the frequency from f2 in the 
high energy limit. Note that A 2 /Q ~ 0(1)- Equation (joTj) is hard to solve because of 
the summation on the right hand side of the equation. We thus look for an approximate 
solution in the leading order of ^ in the large Q limit (D, ^> N). It is clear by examining 
the equation that the last two terms on the right hand side of the equation makes a higher 
order contribution, and can be neglected in a large field expansion. Dividing both sides of 
Eq. (|5Tj) by u> — Ak ~ Ak ' and performing the summation over k', we obtain 



((I k ri + S-;S + )) 



1-P* A 

8^ k , u + 



A 2 



G] (u) 



(52) 



One can solve for ((/ I k ,S ; following a similar approach, 

1 -P 2 



((I-I+S-;S + }} L 



8tt 



-A k ' E - 



A\ 



fe w H 2 — 



(53) 



We then substitute Eqs. (|52|) and ()53|l into Eq. (joT|) to obtain the solution for 
({I k I k ,S~;S + )) UJ . Finally we use Eq. (j4T?j) to find the Green's function in the high en- 
ergy limit. The result is 

1 



1-P 2 



-Si(cj) 



i-p 2 



[E 2 (cD) + S 3 (cD)] 



(54) 



where 



Si(«) = E- J 
it,*;' u + 



A 2 A 2 



Ak-Ay > 



(55) 
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FIG. 5: The self-energies (left panel) and ^2(10) (right panel) as functions of to = to - CI 

- AijCl- In these numerical calculations, we have assumed that the hyperfine coupling constant 
takes the form as in Eq. Both the real part (dashed line) and imaginary part (solid line) of 
the self-energies are plotted. 



A 2 A 2 A 2 
k,k',k" H o — )\ u 2 — ) 



and 



A 2 A 2 A 2 

v (r,\ — \ ^-k^k'^-k" f K7 \ 

| A k -A k , Yr A k „-A h , y ' ^ 

k,k',k" H 2 )\ U ^ 2 ) 

Because £i(<D) has two summations over nuclear spins, £i(u)) ~ A^. For the same reason, 
£2(0;) and £3(0)) are proportional to iV 3 . In the above derivation, we have neglected the 
difference between Ylk"(k,k') an d Hk"i which only gives an order unity contribution to the 
three self-energies. The final expression of G±(co) we have found in Eq. (j54j) is of order 
0(1), in contrast to the solution in the low-energy limit. 

The evaluation of self-energies Si (a;), £2(0)), and £3(0)) is quite complicated. In Appendix 
|D]we first show that £2^) = ^O^)- After converting the summations into integrals, we 
then perform two-dimensional and three-dimensional integrations for £i(u)) and £2(0)). In 
Fig. El we plot the real and imaginary parts for both of the self-energy terms as a function 
of the shifted frequency Co. We can see that there is a cusp at to = for the imaginary parts 
of both £i(o>) and T, 2 {ijj). Let us first take a closer look at £i(o>). Notice that Im£ 1 (u)) is 
negative for the whole region of to. This ensures that the spectral function p(Co) is always 
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FIG. 6: The spectral function p(u>) in the high energy limit as a function of the shifted frequency 
u>. We have plotted the spectral function for 9, = 1.5N, 0, = 2N, Q = 2.5N, and = 4.5JV. In 
each of these cases, we have calculated piCo) with two approximations. In the first one, we have 
evaluated p(u>) with only Si ((D) (dashed line); and in the other one, we have used both Si ((D) and 
S 2 ((D) (solid line). The lower figures clearly show the success of our large field expansion method. 
These plots also illustrated the idea that the original delta function in the spectral function is 
broadened to a continuous spectral function after the virtual nuclear spin flip-flop is included in 
the calculations. 

positive. Like lma n (u) [see Eq. ()42|) and Appendix C], which is nonzero when — ~ < oj < |, 
ImEi(u)) does not vanish when — | < uj < |. The behavior of E 2 ((D) is very different from 
Si ((D). Neither the real nor the imaginary part of S 2 (u)) is symmetrical about Q — 0. On the 
other hand, ImE 2 (d)) changes sign near uj = 0. This leads to no contradiction because E 2 (u)) 
is a higher order correction term to Ei(d)), and E 2 (d)) by itself has no physical meaning. Fig. 
El shows that the absolute values of E 2 (d)) are generally larger than those of Ei(cD), these 
self-energy terms are multiplied by ^2 and g^3, respectively, in Eq. (J54|) . so that Ei((D) is 
still the leading order contribution. 
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So far the two-pair correlation function has been neglected when deriving the EOMs. 
If we keep the two-pair correlation function and all the other corresponding functions of 
the same order, but neglect the three-pair correlation functions, we would find two more 
higher-order self-energy terms which are proportional to p^yr and . In principle, we 
can perform the large field expansion to higher and higher orders, and eventually, we should 
arrive at a finite geometrical series of We would get a finite instead of an infinite series 
because there is a finite number (N) of nuclear spins in the quantum dot, and there can be 
at most iVj_ (iVf) pairs of flip-flopping nuclear spins if N± < (TVf < iVj). This gives rise to 
a natural cut-off of self-energy terms in our large field expansion. However, even extension 
of the calculation to the order of and is already very complicated. Nevertheless, 
the approximation with the first two self-energy terms should be accurate enough as long as 
the total effective magnetic field Q is large compared with N (representing the Overhauser 
field produced by a polarized nuclei spin reservoir). 

The validity of our large field expansion method is illustrated in Fig. where we have 
computed and plotted the spectral function p(u) for various Q. To compare, we have cal- 
culated p(ui) using only the self-energy Si (a;) and both Si(u)) and £2(0)). in the two upper 
panel figures, where Q is not very large, the contribution of S 2 (<I>) is comparable to Si (a)). 
In this case, one needs to calculate more higher order terms of self-energy in the expansion to 
achieve convergence. When Q = 2.5N, the difference is already not so significant; and when 
Q = 4.5iV, we can clearly see that the contribution of S2(u)) is negligible. Therefore, our 
results should be accurate enough as long as Q > 2.5N. For unpolarized nuclei, Q = 2.5N 
corresponds roughly to an external field of 5 Tesla's for GaAs dots.-^ 

We have checked the sum rule for the spectral function, Jf^ p(Q)du = 1, in all the 
numerical calculations. We find that the sum rule is always satisfied within 10~ 3 . Fig. El 
also shows that the original delta function 5 [u — Q + 0(1)] without decoherence has been 
broadened to a continuous spectrum after the electron-mediated flip-flop of nuclear spins 
are explicitly considered. The sum rule from the continuous spectrum clearly indicates that 
there is no contribution from a delta-function to the spectral function, so that the decay of 
G±(t) should be complete. 

Once the spectral function is obtained, the electron spin dephasing time T 2 can be es- 
timated as the half- width of the spectrum at uj = peak, i.e., T 2 = 1/Au). The result is 
shown in Fig. where the dephasing time is plotted as functions of the nuclear polarization 
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FIG. 7: The electron spin T 2 time in logarithm scale as a function of nuclear spin polarization 
P (left panel) and the effective magnetic field O (right panel). T 2 is estimated with the spectral 
function by finding the half-width of the spectral peak. It is easy to see that the coherence time 
can greatly enhanced by increasing the nuclear polarization to 80%. We have assumed A = 92 [ieV 
in these calculations. 

P and the effective magnetic field Q. In all the numerical calculations, we have assumed the 
summed hyperfine coupling constant A = J2k Ak = 92 fieV. Fig. [^indicates that to increase 
T 2 , polarizing the nuclear spins is more effective than increasing the external magnetic field. 
In essence, increasing P leads to a reduction to in the phase space for nuclear spin flip-flop, 
while increasing Q leads to increased energy for the intermediate state and thus reduced 
cross-section for the higher-order process. The T 2 time can be enhanced by almost one order 
of magnitude by increasing P from to 90%, while applying higher magnetic fields extends 
the coherence time by a few times. For unpolarized nuclei with Q = 2.5N, we find T 2 ~ 10 
microsecond, which is similar to the decoherence time caused by the nuclear dipole-dipole 
interaction. This comparison needs to be kept in perspective, however, since our numerical 
results are calculated for an effectively two-dimensional QD with Gaussian type wave func- 
tion, while the previous results for dipolar coupling are obtained for a 3-dimensional QD.— 
It would be interesting to explore how dimensionality changes the T 2 times quantitatively. 
In addition, our results are obtained by assuming nuclear spin I = | throughout this paper, 
though all the isotopes of Ga and As nuclei have spin I = |. Exploration of the dimension 
and / dependence of T 2 would be interesting, but goes beyond the scope of the present 
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FIG. 8: The real part of the Green's function G±_ for different regime of T 2 . G± is evaluated 
by performing the inverse Fourier transform with Eq. (|14|) . We have only plotted the envelope 
function of the correlation function. The actual evolution should be modulated by a fast oscillation 
with frequency fi. 



paper. 

Obviously the T 2 time is proportional to the unit N/A. Therefore a larger quantum dot 
[with a larger number of nuclear spins (iV)] has a longer T 2 time. In the bulk limit the 
coherence time becomes infinitely large because the hyperfme coupling is homogeneous so 
that there is no fluctuation of nuclear field. In the other limit of a smaller quantum dot the 
profile of the electron wavefunction becomes sharper because of the strong confinement. One 
thus also anticipate increased T 2 time since the effective nuclear flip-flop is more difficult to 
realize energetically just like in the case of the dipolar coupling.— These arguments indicate 
that there should be a certain point regarding the QD size where the decoherence effect is 
the most serious. However, in our present calculations we would not be able to reach the 
latter limit because we cannot discuss the influence of the distance among nuclear spins on 
the electron-mediated interaction when converting the sum over nuclear sites into integrals. 

The real time dynamics of the transverse electron spin Green's function G±(t) can be 
obtained by performing the inverse Fourier transform using Eq. f)14jl. In Fig. |S]we show the 
evolution of the envelope of G±(t) (the maximum values during each cycle) for various P and 
Q combinations. The actual evolution of G±(t) oscillates coherently with frequency close 
to Q. In the relative small-magnetic-field (<10T) and low-polarization limit, the amplitude 

28 



of G±(t) decays rapidly and completely. On the other hand, the coherence of the electron 
spin can be maintained for a much longer period if P = 0.9 even without a very large 
Q. We should mention that the 90% nuclear polarization could be realizable in future 
experiments using dynamical nuclear polarization by polarized electronic transpor t 30 ! 31 or 
circularly polarized photon o 29 i 32 i 33 through hyperfme interaction. The current experimental 
realizable nuclear polarization in QDs remains less than 50%.— Our approximation should 
still be accurate with such a large nuclear polarization, because the numbers of up and 
down nuclear spins are still in the same order of magnitude at P = 0.9. If we go to the 
extremely polarized limit, the decoherence effect should become negligible as in the case of 
fully polarized nuclei. For the curve with P = 0.9, we also observe that there is a clear 
revival for G±(t) even after 100 microseconds. 

With our approach we can identify the long time asymptotic behavior of G±(t). Again, 
the contour shown in Fig. H|can be used to calculate the asymptotic integrals of Eq. (fTlj) as 
t — > oo with the method of steepest descent. For simplicity we will only consider the effect 
of To evaluate these integrals we need to find the asymptotic form of the spectral 

function. The calculation is similar to what we have done for the low frequency solution. 
On the contours of C*2 and C3 in Fig. HJ we define Co = T — is. On the contours of C\ and 
C4, we replace Co by ^\ — is. Here s is a new real variable. The long time (t — > 00) behavior 
is determined by the asymptotic form of the spectral function as s -> + . On Ci and C3, 
we find 

24S7 2 

P(° T ~ is ) = n 2 {1 _ P 2 )N 2 ( 1 T 3k) + °(* 3 )- (58) 

Performing the integral along the contours C2 and C3 as t — > 00, we find 

1440 2 1 

On the contours C\ and C 4 , — is) °c s 2 - Therefore, the integral along C\ and C 4 

has a contribution of 4- which is negligible compared to 4 at long time limit. This non- 
exponential long time decay behavior is slower than that of the low-energy solution with 
strong magnetic field, where the form of \ is found. The 4 behavior here is obtained by only 
considering £i(u)). The inclusion of ^2(00) is too complicated. However, we would expect 
power-law decay instead of the decay which appears in the solution for the low-energy 
solution in the absence of external magnetic field. Notice that the expression in Eq. (jBT?|) 

29 



diverges as the nuclear polarization goes toward 1. However, this is insignificant since our 
discussions in this subsection apply to partially and unpolarized nuclei so that P is always 
less than 1. Instead, the low energy solution has the right limit when P = 1. This is because 
both the exact solution of fully polarized nuclei and the low energy solution include only 
the real electron-nuclei flip-flop, while the high energy solution found in this section is due 
completely to the higher-order processes, which are only significant when P < 1. 

IV. DISCUSSION 

The calculations given in this paper can be generalized in several ways. Firstly, although 
we have been dealing with nuclei of spin \ exclusively throughout this paper, our study 
can be easily extended to nuclei with higher spin values. The qualitative behavior of the 
spectral function and the long time decay we have discovered here should not be modified 
in any significant way by considering nuclei with larger spins. Secondly, we have assumed a 
relatively simple form of the hyperfine coupling constant as a function of the position of the 

k 

nucleus at r^, i.e., = e~jv, to simplify the algebra. This form of corresponds to a two- 
dimensional quantum dot with Gaussian electron wavefunction. However, our derivations 
and approximations are not limited by this choice of Ak and are applicable to any form 
of quantum dot with arbitrary electronic wavefunction. In short, all our results before 
converting the summations into integrals are correct for a general electron wavefunction. 
Generalization to a 3-dimensional quantum dot is more complicated mathematically because 
one would encounter integrations that cannot be performed analytically when converting the 
sums over nuclear lattice sites to integrals. 

A more challenging and interesting generalization of the present work is to consider 
smaller effective magnetic fields, for example with f2 ~ 0(v / iV). It is well-known that the 
thermal statistical fluctuating nuclear magnetic field has this order of magnitude.^ In such 
a case, our large N expansion can still be applied, because y/~N is still much larger than A^. 
The Zeeman energy mismatch determines that the direct electron-nuclei relaxation effect is 
still negligible. However, the large field expansion we have developed in this paper loses its 
power in this situation since the expansion parameter ly ^> 1. Basically one has to consider 
more higher-order correlation functions with more pairs of flip-flopping nuclear spins. We 
expect this to be a very difficult task for analytical calculations. 
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We have studied the Non-Markovian dynamics of electron spin with the nuclear spin 
reservoir in an initial product state. It has recently been found that the electron spin 
dynamics shows different behaviors for randomly correlated initial nuclear spin states or 
ensemble averaged initial states.— Our analytical results can be directly used to investigate 
the electron spin decoherence by summing over the product states in the ensemble. Our 
treatment can also be applied to study the Non-Markovian spin dynamics of spin-boson 
model, where a single spin interacts with a reservoir of bosons at thermal equilibrium in the 
initial state. It would also be interesting to extend the approach presented in this paper to 
study the dephasing time of two-electron spin states in double QD, where both the T 2 * and 
T 2 has been measured.— 



V. CONCLUSIONS 



To summarize, we have obtained the decoherence of a single electron spin confined in a 
quantum dot with large effective magnetic fields (including both external field and nuclear 
field) by computing the Green's function G±(t) for the electron spin using the equation- 
of-motion method. We have solved this problem for three types of initial nuclear spin 
states, the fully polarized case, the almost fully polarized case with one nuclear spin flipped, 
and the more general partially polarized and unpolarized cases. We have found that the 
spectral functions have quite different properties at the low frequency (u ~ A/N) and 
the high frequency region (u ~ A). By comparing the exact solution of fully polarized 
nuclei and the solution of almost fully polarized cases, we demonstrate that the virtual 
nuclear spin flip-flop leads to a new feature in the high energy limit. Our studies of electron 
spin decoherence in the low energy region recover previous results, which did not include 
these higher-order processes. Namely, the decay amplitude is of the order O (j^V More 
importantly, we have also obtained the solution of partially polarized and unpolarized nuclei 
by considering the indirect nuclear spin flip-flop explicitly, helped with a large field expansion 
method. We find that the electron spin T 2 time is very sensitive to nuclear polarization when 
P > 0.6, and the electron spin coherence time can be enhanced by 10 times if the nuclear 
polarization is increased to 90%. We also find that this decoherence is complete, and the 
long time asymptotic behavior of the Green's function representing the transverse electron 
spin dynamics is 4. 
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APPENDIX A: EXACT SOLUTION OF FULLY POLARIZED NUCLEI IN 
JORDAN-WIGNER REPRESENTATION 

For spin one-half nuclei, we can introduce the Jordan- Wigner representation:— 

S+ = d, S~ = d + , S z = --n d 

1+ = e -^a k , 1,7 = a+e*"*, F k = \- n k , (Al) 

where = d + d and n k = a^a k . All operators obey the standard fermion anti-commutation 
relations: {d,d + } = 1; {a k) a~l} = 1. This representation preserves all spin commutation 
relations except for those of two nuclear spin operators at two different lattice sites. However 
these commutators are not necessary for fully polarized nuclei, because there should be only 
one flipped spin at any moment during the evolution. In terms of correlation functions, this 
means we won't encounter any correlation function containing two nuclear operators so that 
there is no problem of ordering different nuclear spin operators using the fermion operators. 
Transforming the original Hamiltonian of Eq. (JTJ into the new representation, we arrive 

at 

+ Ey {*U + > (A2) 

k 

where A = J2k A k . We have ignored a constant term in the derivation. The quartic interac- 
tion term J2k A k n k n,i in i?jw is dropped because J2k A k n k and cannot both be 1 (down 
state). The new Hamiltonian i^jw is in the form of the non-interacting Anderson impurity 
model,— describing a localized state (electron spins) interacting with semi-continuous states 
represented by the nuclear spins. This is the key feature that we focus in this paper: a 
single electronic spin interacts with many nuclear spins with different strengths. The non- 
interacting Anderson impurity model can be solved exactly. Specifically, the exact solution 
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of -ffjw for the Green's function {{d; c? + ))^ is^ 

((d; d + )) u = J- p-. (A3) 

+ CJ + 2" - 4 Efc — *r 

The spectral function (— lm((d; d + )) LU /Ti) that represents the overlapping of the localized 
state with the exact eigenstate is the same as what we have found in Section III. 

APPENDIX B: EOMS FOR THE POLARIZED NUCLEAR RESERVOIR WITH 
ONE-FLIPPED NUCLEAR SPIN 

The EOMs are generated by computing the commutators in Eq. The highest order 
correlation functions that survive are those which involve two spin-lowering operators, either 
S~ or I k . All higher order functions with more spins flipped to down direction vanish because 
the total angular (J2k Afc-^I + $ z ) momentum along z direction (of the external field) is a 
constant of motion. Explicitly, the equations are: 



A 3 

w-Q- S (cj) + A k + 



4(^ 2 - AH A) 



(M-;5 + )), = (* |^|*o> 



1 x - A k A k ,V kk ,(uj) _ 1 „ A k A k ,V k , k (u) 

a; - 4f ) = S + )) w + ^V w (w), (B2) 

+ E ^«^r^;5 + »«, (B3) 

k"(k,k') 

-^<&>*;^>>»- E ^<<4^4i;S + >>., (B4) 
z k»{k,k>) z 

+ ^((4^ + 4"'; - ^f-m#&-; S + )) w . (B5) 

Together with Eqs. (fTUJ). these equations form a closed set. The full solution of these 
equations is mathematically intractable because of the number of equations involved. In 
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Section III, we find approximate solutions in the low frequency (u ~ 1) and high frequency 
(a; ~ Q) regions using expansion. 



APPENDIX C: EVALUATION OF a n (to) 



In this appendix the real and imaginary parts of the terms a n (uj) [defined in Eq. 
appearing in the low energy solution of partially polarized and unpolarized nuclei in Section 
III are calculated. The spectral function of the electron spin correlation function calculated 
from cr n (uj) can be used to obtain the renormalized spin precession frequency and decay 
of electron spin coherence. There are three steps in these calculations. First we perform 
analytical continuation by replacing uj with u + i0 + to obtain the retarded expressions. We 



then use the relation 



i 



P- — in5(x) to separate the principle values and the imaginary 
parts. Finally we convert the summation over the nuclear sites into an integral, J2k ~^ Io° dk. 
The validity and efficiency of this conversion have been discussed before.— 
Recall that the summation 02(0;) is 

A k x , A k 



°2M = 



k u - — 



E 



(CI) 



k *~ 1 2 

Using the procedures described above we find that the real and imaginary parts of 02(0;) are 



Re 0-2(0;) = -2N 





1 




1 




In 


1 + — 


+ In 


1 






2uj 


2lu 





(C2) 



and 



Im cr 2 (u;) 



-2Nn < u < i 



2Ntt 



- § < u < 0. 



(C3) 



03(0;) and 04(0;) can be evaluated in a similar manner. We obtain 

2uj + 1 



Re 03(0;) 



-AN + 4Nu\n 



2uo 



(C4) 



Im 03(0;) = — ANttuj, < uj < -, 

2 2 



(C5) 



Re 04(0;) 



-2N - 8Nuj' 



In 


1 


+ In 


1 






1 






+ 2^ 


2lu 





(C6) 
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and 

-8iVvru; 2 < lu < ± 

Im a^co) = { 2 (C7) 

Wtvlo 2 - \ < to < 0, 

Both to = ±| and uj = are branch points for the self-energy. The two branch cuts ([-1/2,0] 
and [0,1/2]) come from different dynamical fields of the electron felt by the nucleus, i.e., 
^ when S z = |, and — ^ when S* z = — |. In contrast only u; = — | appears as a branch 
point in the fully polarized case because S z = 1/2 (rid = 0) makes no contribution to the 
Hamiltonian in IA2I 



APPENDIX D: EVALUATION OF SELF-ENERGY TERMS, Si(u>), E 2 (w), AND 



E 3 (w) 

Using the definitions of the self-energy terms given in Eqs. (|55j). f)56|) and (|57j). and 
converting the summations into integrals, we can obtain analytical expressions of the real 
and imaginary parts of these self-energy terms. The same procedures as we have described 
in Section O are used in the following calculations. We find 



A 2 A 2 

sr f~\ \ " ^k^k' 



k,k' 10 H 2 ~ ' 2 

where we have written A/, and Ay as the integral variables x and The two-dimensional 
integral can be calculated, so that 



f l f 1 xy 
/ dx / dy _ . 
jo Jo uH — ^ 



;di) 



Re Ei (a;) 



2iV 2 
— 3u)ln 



a; + 4u 3 ln 



4cD 2 



1 


+ In 


2ui- 1 








1 4£ 2 




2cD + 1 





(D2) 



and 



Im Ei (a)) 



2iV 2 



4p| 3 - 3|w| + 1 



(D3) 



Again, the imaginary part of Ei(u;) is nonvanishing only when — | < d> < |. Similarly, 

S 2 (o3) = iV 3 /' dx {' dy f 1 dz ——^—— . 

Jo Jo Jo (lu + ^)(uj + ^) 



ReE 2 (c<;) 



AN / ds s 
Jo 



1 + (2lu + s)ln 



, )0 

2cD + s - 1 



(D4) 



2cu + s 



(D5) 
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Im Sofa)) 



I -Hi 

8<I> 3 



2a; | 



-2a) 4 
13 — ■'■ 



f ln|l 



2cD| 



10ii 2 



< a) < \ 



-2a) - 8a)ln|l 
8cu 2 ln ' 



j_ 

2d) 



- 2cD| + if^ln|l 

14£ 2 _|_ 8£ 3 
3 -t " 3 



2a) I 



-In 



4 + 1 , -|<^<o, 



2& 1 ' 2 

Repeating the calculation for S 3 (o)), we find that E 3 (a)) = S 2 (a)). The integral in Eq. 
flD5|) can be computed numerically. However, the calculation is non-trivial because of the 
singularities appearing in the expression. Alternatively, the integral can be done analytically 
using Maple. The result, which is too complicated to be presented here, is a sum of terms 
that include the dilog functions defined as /diiog(^) = Jo hit/ (l—t)dtM Numerical calculation 
of the real part of £ 2 (a>) using the analytical expression obtained with Maple is very accurate 
when we check the sum rule of the spectrum function numerically (see discussion in Section 
III). 
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